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In a recent computational study, we found highly structured ground states for coarse-grained 
polymers adsorbed to ultrathin nanowires in a certain model parameter region. Those tubelike 
configurations show, even at a first glance, exciting morphological similarities to known atomistic 
nanotubes such as single- walled carbon nanotubes. In order to explain those similarities in a sys- 
tematic way, we performed additional detailed and extensive simulations of coarse-grained polymer 
models with various parameter settings. We show this here and explain why standard geometrical 
models for atomistic nanotubes are not suited to interpret the results of those studies. In fact, the 
general structural behavior of polymer nanotubes, as well as specific previous observations, can only 
be explained by applying recently developed polyhedral tube models. 



I. MOTIVATION 



II. REVIEW OF THE TRADITIONAL 
APPROACH TO NANOTUBE GEOMETRY 



The considerations of the importance of exact ge- 
ometric calculations when dealing with curved nanos- 
tructures that we now present arose from observations 
we made during a recent study [l|-[3|. There, we in- 
vestigated polymers adsorbed at ultrathin nanowires by 
means of Monte Carlo simulations applying a common 
coarse-grained bead-stick model. For very high adsorp- 
tion strengths we found, independently of the effective 
radius of the nanowire, well-ordered tubelike monolayer 
ground-state structures for that system. Those polymer 
tubes are formed by aligned helical monomer strands and 
possess different chiralities for different tube radii. Sim- 
ilar structural behavior is known from several tubelike 
atomic structures with applications in nanotechnology 
such as single walled carbon nanotubes 0-0 • 

In order to reveal the morphological similarities be- 
tween the polymer monolayer tubes and atomic nan- 
otubes, we first review the common geometrical view of 
these structures (Sect. [III). We show in Sect. IIIII that 
these approximations are not suitable to explain our find- 
ings and how it should be corrected. Finally, we show 
in Sect. HVl how exact geometrical calculations of curved 
discrete tubes can provide the link between previous re- 
sults [l| and results from additional exhaustive and de- 
tailed computational studies of ground states of polymer 
nanotubes and real- world atomic nanotubes m • 



A common conception about single- walled carbon nan- 
otubes is that they are 'built of rolled up and optionally 
tilted graphene sheets which are carbon atoms crystal- 
lized in a monolayer honeycomb lattice, as depicted in 
Fig. [TJ The standard geometrical description is hence 
based on this corresponding unzipped planar represen- 
tation, which is uniquely defined by a wrapping vector 
Ch pointing from an atomic position to its next periodic 
copy (see Fig. [2]). This vector can be represented as a lin- 
ear combination of two base vectors ai and a2 and two 
integer numbers n and m: Ch = nai + ma2. Conse- 
quently, the vector (n, m) is commonly used to classify 
carbon nanotubes. [5] 



In this common picture, the radius r 



{n,nn) 
classical 



of a carbon 



nanotube is calculated by identifying the length of Ch 
with the perimeter length of the tube: 



2 Trr, 



(n,m) 
classical 



(1) 



where a is the edge length in the lattice or the bond 
length between carbon atoms, respectively. The wrap- 
ping angle ^classical defined to be the angle between Ch 
and ai: 
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Although visualizations sometimes lead to premature 
assumptions, they are obviously and doubtlessly quite 
useful and instructive for the imagination and interpre- 
tation of scientific data, as Figs. [1] and [2] exhibit. Fig- 







FIG. 1. Snapshot sequence from an animation illustrating the unzipping of a carbon nanotube. 





FIG. 2. Definition of the wrapping and the base vectors in three different unzipped planar honeycomb structures. Left: (6,0), 
middle: (6,2), right: (6,6). See text for nomenclature. 




FIG. 3. A triangular nanotube with (n, m) = (3,0), shown 
from different perspectives. See text for nomenclature. 



ure[T] shows snapshots from an animation made with the 
animate-package Q for The input picture se- 

quence was created using the latest Atomistic Simulation 
Visualization software AViz p^, JJj . 



III. THE POLYHEDRAL MODEL AND EFFECT 
OF CORRECTION TERMS 

Comparing tubelike ground states of adsorbed poly- 
mers and carbon nanotubes, one first notes that the 
polymer does not crystallize in a honeycomb lattice, but 
rather in a triangular lattice. See Fig. [3] showing a tri- 
angular (3,0) tube, for example. However, this does not 
change much in the above described picture. As the lat- 
tices are dual in a sense, i.e., one can imagine the sites 
of the triangular lattice residing in the vacancies of the 
honeycomb lattice or centers of the hexagons, one can 
obviously use the same notation and calculation as in- 
troduced above. The only difference is that a has to be 
scaled by a factor of ^/3. 

However when one compares numerical details, one 
notes almost immediately that the above described com- 
mon approximation does not provide a suitable model 



to describe our computational results. To illustrate 
this, we present a comparison of radii in Fig. HJ The 
top row (upside-down solid triangles) shows the solu- 
tions of Eq. ^ for all n and m values corresponding to 

^classical ^ bottom row (upright open triangles) 

represents radii where we find defect-free chiral ground 
states of polymer tubes in our simulations. Obviously, 
there is no apparent link between both sequences. We 
will illustrate (assuming that the bond length between 
monomers or atoms remains fixed) how the deviations 
can be explained using two simple examples. Consider 
first a triangular (3,0) tube as shown above in Fig.[3l i.e., 
a tube where three bonds and three monomers lie in a 
plane perpendicular to the tube center. With a = 1, one 

gets ^cCi'sTicai ~ 3/(27r). However, the three bonds form a 
triangle where the sum of the edge lengths equals 3 and 
the radius of the circumcircle reads r^^a^^ = a/3/3. As 
a further example take the (3,3) triangular tube, where 



IChl = 3^/3 and hence 



^(n,m) 
classical 



3V3/(27r). On the 



other hand the six corresponding bonds form, projected 
to a normal plane, a hexagon with edge length cos(7r/6). 
The radius of the circumcircle equals the edge length in 
a hexagon and therefore r^^act^ = cos(7r/6). In both ex- 
amples, obviously r^Ci'sTiial ^ ^ex^Tt^- 

In fact, when zipping a two-dimensional discrete hon- 
eycomb structure, the geometry changes due to the intro- 
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FIG. 4. Possible radii calculated with usual formula (solid, 
upside-down triangles) and radii of defect-free ground-state 
structures found in simulations (open triangles). 
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duced curvature. Such curved structures should hence be 
described by applying a suitable three-dimensional model 
rather than an effective two-dimensional one. For tri- 
angular tubes and under the assumption that the edge 
lengths remain fixed [12i], such a model is, for example, 
the polyhedral model introduced for the description of 
idealized boron nanotubes [l3|. Calculations within this 
model lead to transcendental equations for both, the ra- 
dius of the nanotubes and the chiral angle, which must be 
solved iteratively. In Fig. [5l we compare the radii of tri- 
angular nanotubes calculated using the classical and the 
polyhedral model as well as the difference for all n and 
m. In Fig.[5fa), where the values of the radii are ordered 
with respect to n and m, those differences do not seem 
to be of importance as the quahtative structure of the 
curve is the same for both models. But if one sorts the 
same data with respect to r, as shown in Fig. [5fb), the 



effect of the deviations becomes apparent. As indicated 
exemplar ily by the vertical dotted lines, there are inter- 
vals where the differences between the radii of two, in 
general completely different, nanotubes is smaller than 
the difference between the radii for the same tube cal- 
culated with the different models. Hence, one can not 
really resolve the hnk between radius and structure of 
nanotubes by using a simple planar ansatz.p^ 

It should be mentioned that when expanding the rela- 
tions for radius and chiral angle [13], the first term of the 
expansion is indeed the same term as one obtains when 
applying the planar model (cp. Eqs. ([T]) and ([2])). How- 
ever the higher order terms are relevant and must not be 
ignored at least in computational studies. [HI 



IV. STRUCTURE OF MONOLAYER 
POLYMER TUBES 
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FIG. 5. Comparison of radii of triangular nanotubes calcu- 
lated in the classical and ind the polyhedral model, (a) Data 
sorted with respect to growing n. Data points between (n, 0) 
and (n + 1,0) correspond to {n,m) with m G (0...n]. (b) 
Data sorted with respect to r. 



If we apply the correct polyhedral model for triangu- 
lar nanotubes, we find indeed a perfect match between 
calculated observables and those found in simulations, in 
contrast to the situation earlier depicted in Fig. |4j We 
show the results of both, calculation and results from 
simulations, in Table HI In the second column, calculated 
radii using the polyhedral model [l3| are given, in the fol- 
lowing three columns simulational details can be found 
and in the last two columns we list the calculated wrap- 
ping angles using a suitable polyhedral model for carbon 
nanotubes [l^ [l7j . 

In order to facilitate the simulations and to obtain 
more precise data, we adapted the model with respect to 
the given problem. In contrast to our recent study [1] we 
introduced flexible nonelastic bonds between monomers 
modeled by the FENE potential and changed the non- 
bonded Lennard- Jones potential such, that the equilib- 
rium distances of both interactions match. Further- 
more we initialized the simulation with a configuration 
where all monomers have the same predetermined per- 
pendicular distance to the nanowire and allowed only 
update moves which not change that distance. In prac- 
tice, we simulated a flexible polymer on cylinder sur- 
faces with more than one hundred different radii (rinput) 
and searched for non-defective ground states, of which 
we measured the chiral angle 0. On the one hand, we 
find the respective polymers tubes indeed for those radii 
calculated from the polyhedral model for triangular nan- 
otubes (cp. columns two and three in Table H]), on the 
other hand we also measure exactly the wrapping angles 
which we calculate for carbon nanotubes (cp. columns 
five and six). 

Hence, the ground states of monolayer polymers nan- 
otubes forming at strongly attractive nanowires can be 
well described by the polyhedral model for idealized 
boron nanotubes and the sequence of wrapping angles 
when changing the radius of the polymer tube is ex- 
actly that calculated for carbon nanotubes. Our simu- 
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polyhedral 


Simulation on 


Cylinder surface ("2D") 


Corresponding 


polymer tube 




output (Ground State) 


carbon nanotube 




' exact 


Tinput 


type 


in ° 


in ° 


type 


(2,1) 


0.51962 


0.477. ..0.532 


3-helix 


20.8. ..17.5 


18.43 


(2,1) 


(3,0) 


0.57735 


0.553. ..0.574 


(3,0) 


0.0 ±0.5 


0.00 


(3,0) 


(2,2) 


0.61237 


0.585. ..0.617 


(2,2) 


31.3. ..29.6 


30.00 


(2,2) 


(3,1) 


0.64526 


0.627. ..0.670 


4-helix 


13.9. ..12.6 


13.57 


(3,1) 


(4,0) 


0.70711 


0.680. ..0.712 


(4,0) 


0.0 ±0.5 


0.00 


(4,0) 


(3,2) 


0.74313 


0.723. ..0.755 


5-helix 


23.8. ..22.7 


23.33 


(3,2) 


(4,1) 


0.78561 


0.765. ..0.808 


5- helix 


11.0. ..10.2 


10.72 


(4,1) 


(5,0) 


0.85065 


0.819. ..0.851 


(5,0) 


0.0 ±0.6 


0.00 


(5,0) 


(3,3) 


0.86603 


0.861 


(3,3) 


29.9 ±0.4 


30.00 


(3,3) 


(4,2) 


0.88462 


0.872. ..0.904 


6- helix 


19.2. ..18.3 


19.01 


(4,2) 


(5,1) 


0.93259 


0.914. ..0.957 


6-helix 


9.0. ..8.4 


8.84 


(5,1) 


(6,0) 


1.00000 


0.967. . . 


(6,0) 


0.0 ±0.5 


0.00 


(6,0) 


(4,3) 


1.00188 


...1.021 


7- helix 


26.2 ±0.7 


25.26 


(4,3) 


(5,2) 


1.03116 


1.031. ..1.052 


7- helix 


15.9. ..15.5 


16.02 


(5,2) 


(6,1) 


1.08319 


1.063. ..1.106 


7- helix 


7.6. ..7.2 


7.52 


(6,1) 


(4,4) 


1.13152 


1.106. ..1.127 


(4,4) 


30.3. ..30.0 


30.00 


(4,4) 


(5,3) 


1.14441 


1.138 


8-helix 


21.7 ±0.4 


21.75 


(5,3) 


(7,0) 


1.15238 


1.148. ..1.169 


(7,0) 


0.0 ±0.6 


0.00 


(7,0) 


(6,2) 


1.18076 


1.169. ..1.201 


8-helix 


13.9. ..13.4 


13.83 


(6,2) 


(7,1) 


1.23600 


1.212. ..1.254 


8-helix 


6.7. ..6.4 


6.54 


(7,1) 


(5,4) 


1.26887 


1.244. . . 1.276 


9-helix 


26.8. ..26.0 


26.32 


(5,4) 


(6,3) 


1.29090 


1.286 


9-helix 


19.0 ±0.5 


19.07 


(6,3) 


(8,0) 


1.30656 


1.297. ..1.318 


(8,0) 


0.0 ±0.6 


0.00 


(8,0) 


(7,2) 


1.33242 


1.318. ..1.361 


9-helix 


12.2. ..11.8 


12.17 


(7,2) 


(8,1) 


1.39027 


1.371. ..1.424 


9-helix 


5.8. ..5.6 


5.79 


(8,1) 



TABLE I. Comparison of calculated observables for nanotubes using polyhedral models and results from computer simulations 
of polymers on a cylinder surface. The rows are ordered with respect to r^^'^^ . 




FIG. 6. A ground- state conformation found in a general study 
of polymers adsorbed to ultrathin nanowires. The structure is 
composed of two competing regions with different chiralities. 

lations are furthermore very precise, we can resolve that 
sequence even for tubes whose radii differ by < 1% (cp. 

^exfct and r^^^l^ or rf^^}^ and r£fjj, which was not pos- 
sible before. 

We can also explain specific results from our general 
study. Figure [6] shows a low-energy conformation with 
tube radius 0.6 < r < 0.7 found for the originally used 
polymer-wire model |1[ . It shows two competing regions 
forming a helix with four strands and a (2,2) structure, 
respectively. For the first one we measure a wrapping an- 
gle ~ 14°. Looking up in Table |T] we find that that part 
corresponds to a (3,1) tube and that this is a direct neigh- 
bor of the (2,2) tube with respect to the possible discrete 
radii for (n,m) tubes. Indeed, it is plausible, that there 
is competition between these structure for radii that do 
not match exactly any r^^'^^ . 



V. SUMMARY 

In this paper we argued that polyhedral models for 
nanotubes are useful for the description of respective 
structures in computational studies. The corrections in- 
troduced by those models compared to the commonly 
used pictures are in general not negligible (see also a 
recent study on the effect of chirality on nanotube vi- 
brations [18|). In particular, the polyhedral model for 
boron nanotubes reflects the findings of monolayer poly- 
mer nanotube structures found earlier. The sequence of 
chiral angles of polymer nanotubes with different radii is 
the same as for carbon nanotubes and provides the link 
between those structures. 
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